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Abstract 

There is evidence for existance of massless Dirac quasi-particles in 
graphene, which satisfy Dirac equation in (1+2) dimensions near the 
so called Dirac points which lie at the corners at the graphene's bril- 
luoin zone. Certain subtle points which are peculiar to odd number 
of space-time dimensions (in this case three) , in the derivation of such 
an equation are clarified. It is shown that parity operator in (1+2) 
dimensions play an intersting role and can be used for defining con- 
served chiral currents [there is no 7 s in (1+2) dimensions] resulting 
from the underlying Lagrangian for massless Dirac quasi-particles in 
graphene which is shown to have chiral Ul(2) x Ur{2) symmetry. Fur- 
ther the quantum field theory (QFT) of Coulomb interaction of 2D 
graphene is developed and applied to vacuum polarization and elec- 
tron self energy and the renormalization of the effective coupling of 
this interaction and Fermi velocity. 

Recently there has been a lot of interest in graphene. On experimental 
side it has some anomalous magnetic and transport properties [1]. On theo- 
retical side its low energy excitations are massless, chiral Dirac fermionsfl]. 
This implies existence of quasi-particles in graphene, which satisfy massless 
Dirac equation in (1+2) dimensions near the so called Dirac points which 
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lie at the corners of graphene's Brilluoin zone. This is indicated by energy 
dispersion around the Dirac point given by [1] 

£ ± ~±^|q|+0(|q| 2 ) (1) 

which reminds one of relativistic energy-momentum relation for massless or 
ultra relativistic particles although Fermi velocity is three hundred times 
smaller than the speed of light. The derivation of such an equation given in 
literature [1,2] is in the form 

icr.Vip = Eip (2) 

where cr = (a 1 , a 2 ) and V = (^fr, gf? ), which has two dimensional form 
of the usual Dirac equation for two component left-handed or right-hand 
massless fermions and not exactly in the form of Dirac equation in (1+2) 
dimensions [see below]. The purpose of this paper is to clarify such points and 
to disscuss the subtleties of odd [in this case three] space-time dimensions 
and develop a quantum field theory of Coulomb interaction of 2D graphene. 

In graphene for the hexagonal layer the unit cell contains two atoms A 
and B, belonging to two sublattices. The lattice vectors are [1] 

ai = |(3,V3), a 2 = |(3,-V3) (3) 

where a ~ 1. 42 A is the carbon-carbon distance. 

The Dirac points K and K have their postion vectors in momentum space: 

In the tight binding approach, the Hamiltonion is given by [1] 

H = -t (4bj + h.c) (5) 

<i,3> 

where aj(aj) are annihilation(creation) operators for an electron in sub-lattice 
A with an equivalent definition for sub-lattice B. £(~ 2.8eV)is the nearest 
hopping energy (hopping between different sub-lattices) 

Let us expand the Hamiltonian (5) at any of two independent lattice 
points[2], which we may take as given in Eq. (3), around the Dirac point K . 
This gives 

nn 1 SH* AB 
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where 

U AB = -t[e ik - ai +^ k - a2 ] (6) 

Writing 

k = K' + q 

and expanding around K' keeping terms linear in |q| , 



This gives 



5H 



SH AB = —[iq x + q y ] (7) 



iq x + q y 

-iq x + q y 

Vfi-q^ + q 2 * 1 } (8) 



where Vf = ^ is the Fermi velocity, a 1 and a 2 are Pauli matrices, q % = 
(g 1 ,^ 2 ), in relativistic notation: q^ = (g ,? 1 ,? 2 ). This does not give cr.q as 
implied in Eq. (2). The Eq. (8) can be put in the form 

5H = Vf^a'q 1 + ia 3 a 2 q 2 } 

= 7°(7-q«/) (9) 

where 7 = a 3 , 7 1 = ia 1 , r j 2 = ia 2 are Dirac matrices in (1+2) dimensions. 
The above is the Dirac Hamiltonian for massless fermion in (1+2) dimensions 
near the Dirac point K . The corresponding Dirac equation is 

Hip = Eip ± , E = ±v f |q| (10) 

or 

= 7°(7HV^ (11) 
Writing do = we can write the Eq (10) in covariant form 

1(7% + + 7 2 W = 

or 

i(7"W = (12) 
which is the Driac equation in (1+2) dimensions for a massless fermion [3]. 
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It is important to remark that if we expand around the Dirac point K, 
we obtain 

SH = u/f-gV-gV] 

= u/fi^VgV -ia 3 a 2 q 2 } (13) 

Now it is known [3] that in 3 space-time dimensions there exists two in- 
equivalent representations for 7-matrices [this is true for any odd number of 
space-time dimensions] : 

7 = a 3 , 7 1 = ia 1 ,^ 2 = ia 2 

7 = a 3 , 7 1 = ia\7 2 = -ia 2 (14) 

We have used the first of these represantations for the expansion around 
the Dirac point K . We take the second representation for the expansion 
around the Dirac point K, whcih is obtained from K by the parity operation, 
defined by the matrix [3] 



A= 1 (15) 




so that detA = — 1. Thus we see that under the parity operation 

q 1 < ► q 1 , q 2 < ► — q 2 and SHk < > 5Hk' 

Taking the two representations mentioned above into account we can 
write the parity conserving Lagrangian[3] as 

£ = ^ + (id)ip + + ip_(id)ip_ 

where 

d = 7% + 7 1 <9i + 7 2 c>2 

d = 7 % + 7 1 5i-7 2 -9 2 (16) 
Parity operation takes the solutions in one representation to the other: 

i>+{x P ) = -7] p 1p_(x) 



where x p = (a; , a; 1 , — x 2 ). It is convenient to transform to new fields [3] 

i>A = ^+ 

4> B = nV- (18) 

The Lagrangian (16) can then be written as [3] 

C = ^ A (i^d^ A + ^ B {i^d^ B (19) 
where under partiy operation now 

^a( xP ) = -iVpl 2 ^B( x ) 
= V p ^b( x ) 

^b( xP ) = - ir lpl 2 ^A(. x ) 

= v p ^a( x ) (20) 

In (1+2) dimensions, there is no 7 5 available as in (1+3) dimensions. But 
still we may define 

^ = tA ^t (21) 

so that the Lagrangian (19) takes the form 

C = V^t^l + ^Ri-fd^R (22) 
This Lagrangian is invariant under two independant chiral transformations 

- e ia ^ L , ^ -> e ia ^ R (23) 

where oll and cur are real, and has thus Ul(2) <8> Ur(2) symmetry. The 
corresponding conserved currents are 

J" R = ^Rl^R (24) 

We can develop quantum field theory for Coulomb interaction of 2D 
graphene in analogy with QED (for another approach see [2,4]). The "free" 
Hamiltonian as implied by Eq. (11) is 

Ho = v f J rf 2 r^ t (r) 7 °7.(-iV)V(r) (25) 

v f J d 2 r^(r)j. {-i V)^(r) (26) 



The instantaneous Coulomb interaction in 2D graphene is 

H I= e - [ dW*2 W l*'f l)w( *'* 2) (27) 
2 J Air |xi - x 2 | 

where n(t, x x ) = ^ x i)V'(^ x i) — ^(^ x i)7°V'(^ x i)- This gives the scat- 
tering matrix element between four fermions 

womb , 2 1 /" > + /" .2„ /" « < /|n(*,xi)n(t,x 2 )|i > 



5 c f f m = -*e 2 - / / rf^x / d 2 x 2 J 1 1 v ' , n (28) 

1 % J J J 47TIXX-X2I 

The integral can be written as J d D x\dP x 2 b{t 2 — 1\) where D = 3 in (1+2) 
dimensions. Expanding ip's into creation and annihilation operators, one 
finally abtains [5] 

iT fi = [u{p' 2 ){~ieY)u{p 2 )- — M / — i u(p[)(-ieY)u(pi) -crossed term] 

Anvf 2[-q 2 + (q.7]) 2 ]2 

(29) 

where ^y, |q| = [— q 2 + (g.r]) 2 ] 1 ^ 2 is Fourier transform of 47r | Xl \_ X2 in two 
space dimensions, i]^ = (1,0,0), so that r\ ^ = 7 and longitudnal photon 
momentum may be taken as = 2(g M — 77. g?/ 4 ), /j, = 0, 1, 2; so that (— I 2 ) 1 / 2 = 
2|q|- 

Thus we may write Feynman rules for 2D coulomb interaction in analogy 
with QED as follows: 

(i) Vertex factor: —ie^, 

(ii) Internal lines 

Photon ("longitudinal") line: ^ — r 

(iii) Spin| massless fermion: ^ 

fa = 7 M p M /i — 0,l, 2, this follows from the Lagrangian given in Eq. (19). 

As an application of these rules, we calculate vacuum polarization II (g 2 ) = 
rj^rj u IP^(g), which arizes from photon self energy due to fermion loop. It is 
given by 

-ill^iq) = (-1)— / -^-rTr[(-ie^)-^(-ieY)-„ , , * , ], 

W 1 J v f J (2n) D lK JJ f + ie y J f+2( J q-q.r ] ) fa + ie h 

(30) 

The integration will be done by using Feynman parametrization and dimen- 
sional regularization. The singular terms combine to give (1 — y)T(l — ^) = 
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r(2-f ), i.e. finite answer, which however cancels with a similar term in the 
trace. The net result is 

n"V) = f -^-w(D + hr(2 - 5) Z 1 ^v^T^)^ (3i) 

^/(2tt)- 2 2 Jo 2 A 

where (D + |) arizes from the trace for odd space time dimensions, A = 
\{q.r]) 2 — q 2 }^ = |q| and 

V = -W((q-Vf - <? 2 ) - W + 8q-vW + 17V - (^ W) (32) 
Using D = 3, we finally get 

H(g 2 ) = = ^-(A) 

e 2 . 7r 

|q| 7T 



47TU/ 

= ff||q|, (33) 

in agreement with the known result [6], where g = is the dimensionless 
[in units fi = 1] effective coupling constant. 

We now discuss how the vacuum polarization renormalizes the interac- 
tion coupling constant. For this purpose we consider the response of charged 
fermion to an externally applied field [Coulomb potential in two space di- 
mensions in momentum space isA Q = — ^ = 'H fl A v , where A v = (— ^y, 0)], 
namely 

u{ier)u^A v (34) 



which is modified to 



U{%ei )M 2[A + II(g 2 )] ^ 
^ n. 1 



2 | q r M iTf 



= (35) 



47TU/ 

"2TqTl + f 

where we have used Eq. (33). Thus the renormalized g, often written as g sc 
[6] is 

9sc = -~Ta* ( 36 ) 

1 2 
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The implications of this result are discussed in [6]. It is instructive to calcu- 
late fermion self energy which is given by 

- ffi = £'(S*H^)i<^>^ (37) 

After making the Wick rotation 1° = il E (E for Euclidean matric) and carry- 
ing out the l E integration (which occurs only in fermion propagator), which 
gives 7T|j|, the rest of integration is in D = 2 dimension. Using Feynman 
parameterization 

» = n I x ~ 1/2{1 ~ xrl/2 ax + b(i-xy (38) 

and the dimensional regularization one obtains 

ie 2 A 2 

-iS(p) = -- (- /b-P fi^T-a ( 39 ) 

LbirVf p 



where we have used 



1 r(l-D/2) /1\ 1_D/2 



( 47r )(i?/2-i) r(x) 



'2 

(InL + 7 — ln47r) 



with L — |p| x(l — x), - signifies the ultraviolet log divergence: In A 
Noting that 



- tri ■ v h = i -v = - t> + i p 

E(p) = £[-^ + 7 »pO]lnA (40) 
4 |P| 

where p° = t> j |p| is the energy. The usual interpretation of this result is 
that coefficient of second term in parenthesis gives the radiative correction 
to energy 

E = v f |p| Z^ 1 (41) 

where 

Z 2 1 = l + 9 -\n-^ (42) 

4 |p| v ' 

The coefficient of fp gives the renormalization of the electric charge, e — > 
Z 2 e,which however is cancelled by the corresponding contribution from the 



8 



vertex part by the use of Ward Identity. The vacuum polarization correction 
is finite and renormalize g to g sc as given in Eq. (36). The result (j4ip is 
interpreted as renormalization of the fermi velocity Vf [2,4,6]. This equation 
is important for the renormalization group analysis of g and Vf [4,6]. 

In summary we have clarified the derivation of Dirac equation for quasi- 
particles in graphene in (1+2) dimensions near the Dirac points K and K. 
The role of two inequvalent representations of Dirac matrices in (1+2) di- 
mensions [a property of odd number of space-time dimensions] visa vis par- 
ity operation is emphasised. It is shown that the underlying Lagrangian 
for quasi-particles in graphene has chiral Ul(2) ® Ur(2) symmetry. Fur- 
ther Feynman rules for QFT of Coulomb interaction of 2D graphite have 
been given and applied to the vacuum polarization and renormalization of 
effective Coulomb interaction constant and electron self energy which has 
important implication in the renormalization group analysis of g and Vf [6]. 

The author would like to thank Ansar Fayyazuddin for introducing him 
to this subject by giving him the reference for the review article [1] . He 
would also like to acknowledge the hospitality of King Fahad University of 
Patroleum and Minerals, Dhahran where a part of the work was completed. 
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